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GAS MOTION IN A LOCAL SUPERSONIC REGION AND CONDITIONS 
OF POTENTIAL-FLOW BREAKDOWN* 


By A. A. Nikolskii and G. I. Taganov 


For a certain Mach, number of the oncoming flow, the local 
velocity first reaches the value of the local velocity of sound 
(M = 1) at some point on the surface of the body located within 
the flow. This Mach number is designated the critical Mach number 
M cr . By increasing the flow velocity, a supersonic local region 

is formed bounded by the body contour and the line of transition 
from subsonic to supersonic velocity. As is shown by observations 
with the Toepler apparatus, at a certain flow Mach number M > M^, 

a shock wave is formed near the body that closes the local super- 
sonic region from behind. The formation of the shock wave is 
associated with the appearance of an additional resistance defined 
as the wave drag. 

In this paper, certain features are described of the flow in 
the local supersonic region, which is bounded by the contour of the 
body and the transition line, and conditions are sought for which 
the potential flow with the local supersonic region becomes impos- 
sible and a shock wave occurs. 

In the first part of the paper, the general properties of the 
potential flow In the local supersonic region, bounded by the con- 
tour of the profile and the transition line, are established. It 
is found that at the transition line, if it is not a line of dis- 
continuity, the law of monotonic variation of the angle of inclina- 
tion of the velocity vector holds (monotonic law) . An approxima- 
tion is given for the change in velocity at the contour of the body. 
The flow about a contour having a straight part is studied. 

In the second part of the paper, an approximation is given of 
the magnitudes of the accelerations at the interior points of the 
supersonic region. With the aid of these approximations, it is 
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shown that for profiles convex to the flow the breakdown of the 
potential flow, associated with an increase of the Mach number of 
the oncoming flow, cannot be due to the formation of an envelope 
of the characteristics within the supersonic region. 

On the basis of the monotonic law, the transitional Mach 
number M is found, beyond which the potential flow with local 
supersonic region becomes impossible. 


I - PROPERTIES OP POTENTIAL FLOW WITH LOCAL SUPERSONIC REGION 
1. Properties of Supersonic Flow Bounded by Solid Wall 

and Transition Line 

The flow in the supersonic region about a wing profile of 
usual form with a flow about it of a Mach number higher than 
critical is considered. In accordance with these conditions, two 
types of such flow may occur: 

(a) Transition from supersonic flow to subsonic flow. - In 
this case, the supersonic region is bounded by part of the solid 
contour AB and the transition line A = 1, which Is indicated 
by a dotted line, as shown in figure 1. The magnitude A => w/a#, 
where w is the absolute magnitude of the velocity and a* is 
the critical velocity (the velocity at which the flow velocity is 
equal to that of sound) ; 

(b) Transition from supersonic flow to subsonic flow with 

aid of shock wave. - In this case, the supersonic region is bounded 
by the solid contour AB, the transition line A = 1, and the 
shock wave CB (fig. 2). 

Through each point of the supersonic region, two character- 
istics of different families pass,. The minimum angle between the 
direction of the velocity at a given point and each of the char- 
acteristics is equal to the Mach angle a. 

If the velocity vector is rotated by the angle a in a 
counterclockwise direction, the direction of the vector will coin- 
cide with one of the characteristics that shall be denoted as the 
characteristic of the first family in contradistinction to the 
other characteristic, which shall be denoted as the characteristic 
of the second family. 
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In case (a), each of the characteristics drawn from, an arbi- 
trary point P of the contour has a point in common with the 
transition line X = 1 (fig. 1). 

In case (b) , all the characteristics of the second family end 
on the transition line, but not all characteristics of the first 
family possess this property. 

By displacing the point P on the contour in the direction 
toward the base of the shock wave B, it is found that originating 
from a certain point D of the contour, the characteristics of 
the first family no longer end on the transition line but intersect 
the shock wave (fig. 2). The essential difference of the flow with 
local supersonic region from other mixed flows (for example, flow 
from a Laval nozzle) lies in the fact that from each point of the 
contour there originates at least one characteristic that ends on 
the transition line. 

The region ADCA, in which both characteristics drawn from the 
points of the contour end on the transition line. Is denoted as 
region Ij the region DCBD, in which only the characteristics of 
the second family fall on the transition line, is denoted as 
region II. 

In case (a), the entire supersonic region coincides with 
region I, hence, all theorems derived in this paper and the 
approximations for the flow in region I, obtained from the single 
assumption of the ending of the characteristics of both families 
on the transition line, will also hold for case (a). 

Hereinafter, there will often be considered, together with a 
certain point P of the contour, simultaneous points on the 
line X a 1, which are the ends of the characteristics originat- 
ing from point P. The end of the characteristic of the first 
family yill always be designated by the same letter as the point 
of the contour but with subscript 1 and by the sign *. Similarly, 
the end of the second characteristic Will be designated by sub- 
script 2 with sign *, as shown in figures 1 and 2. 

In the following derivations, the fundamental magnitudes are 
the inclinations of the tangents at an arbitrary point P of the 
contour to the axis of abscissas, denoted always by 0^;' the 

inclinations of the velocity vector on the tr ans ition line at 
points P 1# and P 2# are denoted by and 0 2 *, respectively. 

All three magnitudes 0^, 0-^,, and 0 2# are functions of the 

arc length of the contour. 



4 


KACA TM No. 1213 


* 


An arbitrary point P on the contour in the supersonic region 
is considered. The characteristic of the first family PP2* is 
represented in the plane of the hodograph of the velocity by an arc 
of an epicycloid of the second family P'Pg*' (fig* 3). Let 

0 = - & k be the polar angle between the points Pg#' and 

p l* ' * For "kke magnitude X at the point P of the profile, the 
following equation is obtained: 

\ = ffeg* - e k ) (1.1) 


where X = f (0) is the equation of the epicycloid so set up that 
the equation f(0) = 1 holds. This equation, as is known, has the 
form (reference 2) 


6 = - \ arc sin & " (X -1 ) “ 


arc sin 


(x + D ^ - xj + - 1) I 


( 1 . 2 ) 


In determining 0 2# from the given X and 6^ 


0 2* * e k + <p00 


(1.3) 


where cp is the function reciprocal to f. The graphs of the 
functions f (0) and cp(X) are given in figure 4. 


By considering the characteristic of the first family starting 
from point P (fig. 3), 

X = ffSfc - = 0 k “ ®P( X ) (1.4) 


* 


where f and 9 are the same functions as in equations (1.1) and 
(1.3). 


For the points of region I where the characteristics of both 
families end on the transition line, the following equation is 
obtained from equations (1.3) and (1.4): 


« 0 1* + S 2* 
9 k g 


(1.5) 


r 
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Thus 

TEEOBEM 1. - The Inclination of the contour at any point P 
of zone I is the arithmetic mean of the Inclinations of the veloc- 
ity vectors at the points of the transition line, which are the 
ends of the characteristics starting out from point P. 

Equation (1.3) and the second equation (1.4) permit, for a 
certain distribution of the velocities on the given contour, deter- 
mination of the inclinations of the velocity vectors at the ends 
of the characteristics lying on the transition line. 


2. Monotonic Law of Change in Angle of Inclination 
of Velocity Vector on Transition Line 


If in the equations of an adiabatic gas, the nondimens ional 
velocity X and the angle of inclination 9 of the velocity 
vector to the x axis are taken as the unknown functions, the 
following equations are obtained: 



*( 8ln 



= o 


sin e|^ - 
ox 


cos 9^^ + X cos 0“ + X sin 9^^ — 0 
oy dx oy 


( 1 . 6 ) 


The character of the change of the magnitude 9 along the 
line X <= Xq ■ constant. On one side of this line, let X < Xq 

and on the other side let X > Xq (fig. 5). The normal to the 
line X = Xq is drawn in the direction of decreasing velocity. 

By considering a certain point of this line and by taking the 
direction of the y axis to agree with the direction of the 
normal to this point, from equation (1.6) 

(1 - M^ 2 ) sin |* - Aj Bin e x H + *! 808 S 1 H - 0 


COS 9q ~ - Xq COS 9q ~ - Xq Sin 9i ~ = 0 (1.7) 

f 

where Mq is the Mach number for X = Xq, and 9q is the angle 
between the velocity vector and the positive direction of the 


\ 
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tangent to the line A = Aq at a fixed point of this lipe. By 
eliminating from equations (l.7) the magnitude b9/ba and by 
determining from the obtained equation bS/bs 

be 1 - Ml 2 sin 2 e x 

as “ A x aS U 

In the subsonic flow, 1 - M^ 2 sin 2 > 0 for any 9 and 
because bb/bn < 0, then b9/bs < 0. Thus 9, in this case, 
changes monotonically along the line A = constant. 


By considering the transition line where A^ = 1 and by 

starting from the assumption that the transition line is not a 
line of discontinuity (reference 1), that is, that all the deriva- 
tives bo/bs f b9/bn, b^/be, and bX/bn are f ini te on the tran- 
sition line, the following equation is obtained from equation (1.8) 

= cos 2 9 |^. — (1.9) 

ds on 

and therefore the condition b9/be < 0. Thus 

THEOREM 2. - If a point moves along the transition line so 
that the region of subsonic velocity lies to the left, the velocity 
vector will monotonically turn in the clockwise direction. 

The condition b9/ba = 0 along the transition line, in the 
case where a transition from the subsonic to the supersonic veloc- 
ity occurs, was previously obtained by S. A. Christ ianovich. 

The property of monotonicity is not, however, characteristic 
for supersonic flow. In this case, the inequality 1 - M 2 sin 2 9 > 0 

and the inequality 1 - M 2 sin 2 9 < 0 could hold and therefore the 
value of 5©/ds can change sign. 

From further discussion, it will be clear that the fact 
expressed by theorem 2 determines, to a considerable extent, the 
character of the flow in the local supersonic region and also the 
possibility or impossibility -of the simultaneous existence of the 
subsonic and supersonic flows without change in the potential char- 
acter of the flow. This fact will hereinafter be termed the "law 
of monotonicity on the transition line," or simply the "monotonic 
law." 
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3. Property of Monotonic Change of Velocity Magnitude 
and Its Inclination along Characteristics 
By mak ing use of the results of section 2, it will he proved 

that: 


THEOBEM 3. - If in the supersonic region there is given a 
section of the characteristic of one family, such that the char- 
acteristics of the other family originating from the points of 
this segment end on the transition line , the angle of inclination 
of the velocity vector and the magnitude of the velocity are mono- 
ton ic functions along the given segment of the characteristic . 

By considering a certain point C on the segment AB of the 
characteristic, draw from the points A, B, and C to the line of 
transition the characteristics of the family different from the 
family to which the characteristic AB belongs (fig. 6). The ends 
of these characteristics lying on the transition line are A*, B*, 

and C # , respectively. 

In the plane of the hodograph, points A, B, and C are 
represented by the points A', B*, and C‘, respectively, lying 

on a certain epicycloid y, and the points A*., B # , and C # 

are represented by the points A*’ , B* ' , and C^’ , respectively, 
lying on the circle X » 1. Each pair of points. A* and A*' , 

B’ and B*', and C* and lie on one of the epicycloids of 

the family different from the one to which the epicycloid 7 
belongs (fig. 7). 

As the point C* moves along the epicycloid 7 and its cor- 
responding point Cy_ ’ moves along the circle X = 1, the polar 

angles 8 at the points C' and C*' either simultaneously 
decrease or simultaneously increase because each straight line, 

8 = constant, intersects the epicycloid only at one point. When 
point C moves in the same direction from point A to point B 
along the segment of the characteristic AB, the point C* moves 
along the transition line from point A* to point B* likewise 
in one direction as the characteristics of one family do not 
intersect. In accordance with theorem 2 of the preceding section, 
it follows that the point C* ’ likewise moves along the circle 
X = 1 in tiie same direction. Bence, the polar angle 8 for the 
point C* varies monotonically and, as follows from the properties 
of the epicycloid, the magnitude of the velocity X likewise 
varies monotonically for point C’ . 
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In the supersonic region I about the profile determined in 
section I, both families of characteristics satisfy the conditions 
of theorem 3 so that in moving along any characteristic, the incline- * 

tion of the velocity vector 0 and the magnitude of the velocity X 
vary monotonically. For motion directed toward the transition line 
along the characteristics of the first family, both magnitudes 0 
and X monotonically decrease; whereas the characteristics of the 
second family (the left one), 0 monotonically increases and X 
monotonically decreases. From the obtained property of the mono- 
tonic variation of the magnitude of the velocity and the angle of 
inclination of the velocity vector along each of the character- 
istics of region I, it follows that this region is represented as 
a single sheet on the corresponding region of the hodograph. 

In the supersonic region n, only the characteristics of the 
first family satisfy the conditions of theorem 3. On moving along 
the characteristics away from the. profile, both magnitudes 0 and 
X monotonically decrease. On moving along the characteristics of 
the second family in region II, the property of monotonicity of the 
change in the magnitudes 0 and X does not, in general, hold. 

By making use of the known properties of the characteristics 
in the plane of the flow and in the plane of the velocity hodo- j 

graph, the directions of concavity, as shown in figure 8, are 
obtained in the region of applicability of theorem 3 for small 
Mach numbers. This direction of concavity of the characteristics, 
as follows from the results of S. A. Christ ianovich (reference l), 

takes place only in the case M < =2/ «/3 - x = 1.565. By mak- 

ing use of the results mentioned, it is found that if M = for 
a certain point K in the region of applicability of theorem 3, 
the characteristics have the appearance represented in figure 9. 

Point K is always a point of inflection of the characteristics. 


4. Variation of Velocity at Profile in Supersonic Begion 


It is assumed, as Is usually the case, that in the supersonic 
region the profile is everywhere directed convex to the flow. In 

the supersonic region at the profile, an arbitrary point K and 
the characteristic of the second family KKja* are assumed to 
originate at this point (fig. 10). For point K from equation (l.l), 

X = f (02* - 0^) where f is the fraction introduced according to 

equation (l.l). 


dX 


a(- e k ) 


- - 4 * * * 8 k> 


d0, 


1 + 


2 * 


*1= flj 


* 




( 1 . 10 ) 
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In accordance with theorem 2, ® k ) — 0$ hence, by 

making use of equation (1.10) 

- « k ) (1-ID 

But 0 2 * “ 0 k = cp(?0> where cp is the reciprocal function 
of f. Hence 

- e k ) = f’tpft)] = 

where the differentiation is effected along the epicycloid in the 
plane of the hodograph (fig. 11) . Two points A and B are 
assumed to be on the epicycloid, the polar angle between which is 
equal to d0. The value of X is also assumed to be greater at 
point B than at point A. 

From point A, an arc of the circle X = constant is drawn to 
the intersection with the straight line OB at point C, the 
length of the arc being taken equal to dCT. The angle between the 
tangents to the arc of the circle AC and the arc of the epicy- 
cloid AB at point A is equal to the Mach angle a. Hence, 
dX = tg a da and because da = Xd0, then = Xtg a. 

The graph of the function f ’ [cp(X)] = Xtg a is given in 
figure 5. The inequality (l.ll) therefore assumes the form 

arr^T- (1 - X2) 


As is shown by the preceding inequality, this value depends 
.only on the magnitude X at the given point. By comparing equa- 
tion (1.10) and inequality (l.ll), it is found that inequality (1.12) 
may be expressed as an equation if and only if at a given point of 
the profile dSg*/^®^ =0. If this relation holds true on a cer- 
tain segment AB of the profile (fig. 12), the condition 0 2# = 


constant is obtained. All the characteristics of the second fam- 
ily starting from AB are then represented in the plane of the 
hodograph by the same epicycloid of the second family y (fig. 13) 
because at their ends, lying on the line X = 1, the velocity is 

Ag # (fig. 12) is also 


constant. Hence, the entire region ABB^ 


represented in the plane of the hodograph by the segment of the 
epicycloid A’B 1 , upon which the points of the segment AB will 
also lie. 
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The characteristics of the first family in the region ABB„ A 

2 * 

must, in this case, he straight lines with constant velocity on 
each; that is, in this region a certain Meyer flow (rarefaction) 
must occur. However, as will he shown in section 6, the realiza- 
tion of such flow under the conditions of the problem is impossible. 
Hence, in the relation (1.12), the equality can hold only at certain 
points of the contour. 

In the supersonic region I, the characteristics of the first 
family likewise end on the line A = 1. By applying considerations 
analogous to those previously mentioned, the second inequality along 
the profile in region I is obtained.- 

a(- d a k ) ' - ~ . f 1 - 15 ) 

In the case where the equality in relation (1.13) is attained 
on a certain segment of the profile AB in the region bounded by 
the characteristics of the first family originating at points A 
and B, the transition line, and the profile, a Meyer flow (compres- 
sion) takes place, the straight lines being the characteristics of 
the second family (fig. 14). As will be shown in section 6, such 
flow likewise cannot be realized under the conditions of the problem 
and therefore, in relation (1.13), the equality can be attained only 
at certain points of the contour. By combining relations (1.12) 
and (1.13), the following inequality is obtained on the contour of 
region Is 


dA 


< Xtg a 


(1.14) 


Inequality (1.13) is of little Interest in the case of a shock 
wave that closes the supersonic region because in this case in 
region I, the velocity generally does not decrease so that the rela- 
tion refers mainly to the case of the flow with supersonic region 
without a shock wave. In this case, region I coincides with the 
entire supersonic region. 


The equation for the change in the velocity can also be 
obtained for the case where a certain segment of the profile is 
concave to the flow. Thus, for any point C on the segment AB 
(fig. 15) 


dA 





X = f (®2* “ ®k) 
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But f'(02* " = A tg a, hence 

-jjj 1 < - X tg a (1.15) 

k 

By making use of the law of monotonicity, dj9g # /de k < 0 and 
therefore from inequality (1.15) 

^ < - X tg a (1.16) 

k 

[NACA Comment: equation (1.16) is apparently in error.] 

Thus, on the segment of the profile having a concavity facing 
the flow, the velocity drops and, as follows from relation (1.16), 
its drop cannot he too slow. Inequality (1.16) differs essentially 
from relations (1.12) and (1.13), which show that on the convex 
segment of the profile, the velocity cannot vary too rapidly. 

Strictly speaking, these conditions are true up to those 
changes in the flow that are brought about by an oblique shock 
wave (fig. 15). By considering the characteristic of the second 
family and i ‘ fcs prolongation C^C, it is found that the 

transition through the shock wave in the flow plane corresponds 
to the displacement along the segment of the strophoid C'-^C^ 

in the plane of the velocity hodograph (fig. 16). Hence, the 
transforms of the segments Cq C 2* an<i c i c of 'the character- 
istic lie on two different epicycloids of the second family. How- 
ever, by making use of the fact that at the point C'-j_ the stro- 
phoid and epicycloid have a common tangent and the same radius of 
curvature, it is found that with an accuracy up to small magni- 
tudes of the third order, relative to the changes in velocity in 
the shock wave, these two epicycloids may be considered as coin- 
ciding. Inequality (1.16) holds with the same degree of accuracy. 


5. Flow in Supersonic Eegion Arising from Presence 
of Straight Segment on Profile 

The straight segment AB is assumed to be on the profile In 
the supersonic region. The characteristic of the second family 
CC2* originates from the arbitrary point C, which is assumed to 

be on the straight segment AB (fig. 17). At point C, the 
equation is X* f(0g* ~ ®k) • 
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As point C moves along segment AB in the flow direction, 
the magnitude 0^ remains constant because of the rectilinearity 

of the segment and the magnitude @2*> accordance with theo- 
rem 2, does not increase. Hence, using the fact that f is a 
monotonically increasing function of its argument: 

THEOREM 4. - On a straight segment of a profile in a super- 
sonic region, the velocity in the flow direction does not Increase. 

It is assumed that the straight segment is located in the 
supersonic region I (fig. 18). 

From the ends of the straight segment to the transition line, 
the characteristics of the first family AA]_* and BB^ and of 

the second family AA£ # and BBjj#. are drawn. Let 0^ = 0 q on 

AB. From point C near point A, the characteristic of the 
second family is drawn to the intersection with the character- 
istic AA 1# at point C’. From theorem 3, it is found that 

0(A) >0(C') >0(C), hut 0(A) » 0(C) * 0 O , hence 0(C) = 0 O 
and, as again follows from theorem 3, 0 = Qq over the entire 
segment of the characteristic CC 1 . 

By moving point C from point A to point B, it is found 
that the characteristics AA^* and BBg* necessarily intersect 
and in the triangle ABE, formed by the segment AB and these 
characteristics, 0 = 0 O = constant. Hence, it follows that 
A * constant and the characteristics of both families are 
rectilinear. 

The entire region ABBg^Ag^ is represented in the plane of 
the velocity hodograph by a single epicycloid of the second family, 
as all the characteristics of the second family are represented by 
a single epicycloid passing through the point 0 = 0 q, A « A(A) . 

Hence, over the entire region considered, a certain Meyer flow 
(rarefaction) will take place with straight characteristics of the 
first family. Similarly, the entire region ABB^A^ is repre- 
sented in the plane of the hodograph by a single epicycloid of the 
first family, over which a Meyer flow (compression) will take 
place with straight characteristics of the second family. 

As was previously stated and as will be shown in section 6, 
the realization of Meyer flows under the conditions of the problem 
is Impossible. 
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6. Impossibility of Meyer Flow in Local Supersonic Region 

In the plane x,y, a Meyer flow is assumed to be between the 
characteristics C^ and Cg of the second family with straight 

characteristics of the first family (fig. 19). 

From a certain point A of characteristic Cq, the straight 
characteristic of the second family is drawn to its intersection 
at point B with the characteristic Cg. On the straight line AB, 

the Mach angle a = constant, hence, it intersects the character- 
istics Ci and Cg at the same angle y = it - 2a. 

Set AB «* Z and denote by 3 the angle formed by the straight 
line AB with the axis of abscissas. The straight characteris- 
tic A-|_B]_, infinitely near AB, is considered. The segments of 

the straight lines AB and AiBi are denoted, respectively, by 
dsi and dsg on the characteristics Ci and Cg and the angle 
between AB and AiBi is denoted by d3. Thus, the following 
equations are obtained: 


ds 


- rdg 
1 sin 7 


as 2 = 


(r + l ) dp 
sin 7 


(1.17) 


where r is the distance between point A and the point of inter- 
section of the straight lines AB and AqBq. 

Set the length of the segment AiBi equal to Z + dZ . From 
figure 18, the equation dZ = (dsg - dsi ) cos y is obtained, or 
by making use of equations (1.17) 

dZ = Z ctg 7 dp (1.18) 


In the plane of the hodograph, the characteristic Cq Is 
represented by an epicycloid. 

The directions AB and AqBq coincide with the direction of 
the normals to this epicycloid drawn at the points A' and A* q, 
corresponding to the points A and Aq (fig. 20). The angle 
between the radius vector and the direction of this normal is equal 
to the Mach angle a. Hence, p » a + 6 and dp => da + d8. 

Thus, by making use of equation (1.18) 


~Y = ctg 7 (da + dfi) = - ctg 2a da + ctg y d© 


(1.19) 
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or "by integrating 



In 


A 


'sin 2ao 
sin 2a 


ctg 7 d 9 


( 1 . 20 ) 


where Zq, ccq, and characterize a certain initial point. 

For a > jr/4 and 6 > 9 q , both terms on the right of equa- 
tion (1.20) are positive so that the following inequalities hold: 

l>l0 (i. 21) 

Z Q V sin 2a u */ sin 2a 

If - A -> 1, then a -> jt/2 and Z -»». 

An entirely analogous result is obtained for the Meyer flow 
between the characteristics of the first family with straight 
characteristics of the second family. Thus 

THEOREM 5. - A Meyer flow, between two characteristics of one 
family with straight characteristics of the second family, cannot 
be entirely realized up to the line A = 1 in a finite region. 


If on the characteristic C-j_, represented in figure 19, the 

equality A = 1 is attained at a certain point D, then by making 
use of the preceding result, it is found that the tangent at 
point D to the characteristic C-. is asymptotic for any char- 

A 1 
acteristic C 2 placed relatively to Cj_, as shown in figure 19 . 

From theorem 5, it follows that along the straight part of the 
profile in the local supersonic region, the magnitude of the veloc- 
ity cannot be constant. If such were the case, then between the 
characteristics of the second family originating at the ends of 
this segment a Meyer flow would take place. This realization, 
because of the finiteness of the local supersonic region, is 
impossible as follows from theorem 5. Thus theorem 4 may be more 
definitely stated as follows: 


^-According to S. A. Christianovich, if a Meyer flow originates 
from a straight characteristic on which A = 1, the character- 
istics of the other family do not originate from this straight 
line. 
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THEOREM 6. - On the straight segment of a profile in the super- 
sonic region, the velocity decreases. 


II - CONDITIONS FOE BREAKDOWN OF POTENTIAL FLOW 
WITH LOCAL SUPERSONIC REGION 

The occurrence of a wave resistance for a flow with, large sub- 
sonic velocity about a body is assumed, to be connected, with the 
instant when at any point on the surface of the body a velocity is 
reached, equal to the local velocity of sound. 

■The study of the results of a number of tests has shown, how- 
ever, that the occurrence of a shock wave and therefore the arising 
of a wave resistance sometimes takes place beyond M^. Photographs 
obtained by the Toepler method show in certain cases a local super- 
sonic region extended deep into the flow without any shock wave. 

The necessity for the breakdown of the potential character of 
the flow with local supersonic region, that is, the necessity for 
the appearance of a shock wave, has not yet been theoretically 
demonstrated for a single case. It will therefore be of interest 
to describe certain conditions that are known to be accompanied by 
a breakdown of the potential flow. 


1. Deformation of Contour 

How the deformation of the contour can lead to the breakdown 
of a previously existing potential flow is considered. 

A flow with local supersonic region about a curvilinear con- 
tour for a Mach number of the oncoming flow such that the flow is 
potential (fig. 21) is considered. 

From points A and B of the contour, the characteristics 
of both families are drawn both in the flow plane and in the plane 
of the hodograph. Because of satisfying the condition of monotonic- 
ity of change of the angle on the transition line, the transforma- 
tions of the points A g*, B £*, A]_ # , and B]_* in the hodograph 

plane axe located on the circle A = 1 in the same order as on the 
transition line in the flow plane. The contour is deformed in such 
a way that the segment of the arc AB transforms into a straight 
segment (fig. 22). 
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By virtue of theorem 6, section I, the velocity along the 
straight segment AB can only decrease. By drawing the character- 
istics in the hodograph plane and by making use of theorem 1, it 
is found that the transforms of the points A 2 *> Bg*, A^*, and 
are not located on the circle X = 1 in the order in which 
these points are arranged on the transition line in the flow plane. 
Consequently, the deformation of the contour breaks up the mono- 
tonicity on the transition line. Thus 

THEOREM 7. - The characteristics of the first family that 
originate from the points of a straight segment of the contour can 
never end on the transition line but must fall on a shock wave. 

If the contour deformation considered is effected continuously, 
it is evident that the breakdown of the potential flow occurs before 
the deformed segment becomes straight. This fact permits the con- 
clusion that at a given Mach number of the oncoming flow, any pro- 
file convex to the flow can be so deformed that a new profile is 
again obtained convex to the flow, the flow about which for the 
same Mach number is no longer possible without a shock wave. 


2. Impossibility of Formation of Line of Discontinuity 
Within Local Supersonic Region 

Investigation will reveal whether the impossibility of realiza- 
tion of a continuous flow within the supersonic region for given 
boundary conditions, as expressed in the intersection of the char- 
acteristics of one family and the formation of a supersonic shock 
wave within the region, can be a cause for the breakdown of the 
potential character of the gas motion. 

Karman (reference 3) makeB the assumption that the probable 
cause of the appearance of shock waves is the formation of an 
envelope of the Mach lines in the supersonic zone for a certain 
Mach number of the oncoming flow. 

The envelope of the characteristics of one family, as shown 
by S. A. Christianovich (reference l), coincides with the line 
of discontinuity that is determined as the geometric locus of the 
points at which at least one of the derivatives 3s/3x, 30/By, 

h'X/Sx, or 3>/dy becomes infinite. Hence, to investigate the 
problem of the possibility of the formation of supersonic shock 
waves within the supersonic region on increasing the Mach number 
of the oncoming flow, it is necessary to obtain the value of the 
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acceleration at the points of the supersonic region. It will he 
shown that on the contour of the profile infinite accelerations 
cannot occur. 


Given a certain flow in the supersonic region with continuous 
transition from supersonic to subsonic velocities (sec. I, pt. 1, 
case (a)); assume a certain point K on the contour in the super- 
sonic region and take the x axis in the flow direction (fig. 23). 
Equations (1.6) then become 


(xZ , \ ,, 3© ~ 3X , 3© n 

Thus 

3X d \ 3© _ v dX M 1 dX 
X * 5^ 5 - ■ * ■ ■ I aej 



( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


where k is the curvature of the contour at point K and B is 
the radius of the curvature at this point. By making use of 
inequality (1.14) and equations (2.2) and (2.3) 


3X 

< 1 -V . 

3© 

c 5 

- jj A tg a 

SE 


1 

E 


From equations (1.20) 


3© 


3X 

Sj 

5 g ctg a 

3y 


h 

E 


(2.4) 


(2.5) 


Eolations (2.4) and (2.5) show that the magnitudes 3©/3x, 
3©/3y, 3>/3x, and 3>/3y are always finite on the contour within 

the supersonic region provided that everywhere E / 0. The case 
E * 0 (infinite curvature) at some point is a special case and never 
occurs on real profiles. Therefore, assume that on the contours 
considered, at all points E > Bq where Bq is a certain constant. 

The maximum possible value of the acceleration at the interior 
points of a local supersonic region is computed. The differential 
equations of the characteristics of the first and second families 
in the plane x,y are, respectively. 
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dy = tg (0 + a) dx dy = tg(0 - a) dx (2.6) 


Moreover, along the characteristics of the first and second 
families, respectively, the relations hold (see for example 
reference l) 


O + 9 = 2| = constant O - 0 *» 2tj * constant (2.7) 


vhere 


l -/ Vr~ 


X 2 - 1 


dX 


C(x - D/(x+ D3 x 2 x 


^ ( 2 . 8 ) 


From relations (2.6) and (2.7) 


If - + «> |f * *8(9 - ») 


dx 

51 


(2.9) 


By differentiating the first relation with respect to £, the 
second with respect to rj, and by subtracting one from the other, 
the following equation is obtained for the value of x: 


1 3(0- a) 3x \ 3 2 x 

- V 1 5! + te(e - a) SIS? 


A 


1 S(0+a), 3x + t g(0+a) 

cos 2 (0+a) 3 1 3q 3|3 t] 


( 2 . 10 ) 


From relations (2.7) 

<7 = £ + q 0 «= I - q (2.11) 


By m a k i n g use of these equations and choosing the coordinate 
system so that the direction of the x axis coincides with the 
direction of the velocity at some point A (fig. 24), the relation 
between the partial derivatives of x at the point considered is 
obtained from equation (2.10) 


3 2 x + 1 a* (g) 3x 1 + a* (0) 3x 

3|3£ sin 2a ^1 sin 2a Sri" 


0 


( 2 . 12 ) 
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Let d Z and ds be the lengths of the elementary segments of 
the characteristics of the first and second families, respectively. 
In passing from point A to neighboring points in the directions 
indicated in figure 24, <32 and ds vill he taken as positive 


dx 5x dl 

ag = az ag 


cos a 


dZ 


ax _ ax ds 

<3s dtj 


cos a 


ds 

dtj 


(2.13) 


where the derivative dl/dg is taken along a characteristic of the 
second family and the derivative ds/drj is taken along a char- 
acteristic of the first family. By making use of equations (2.13), 
relation (2.12) can be reduced to one of the two following forms: 


a 


cos a 



. 1 + a '(a) 

sin 2a 






1 + a *(g) 
sin 2a 



1 + a* (a) ds 
2 sin a dt] 

(2.14) 

L4,a'.(pl « . 0 

2 sin a df 

(2.15) 


Along the characteristic g «= constant, relation (2.14) may 
be considered as an ordinary differential equation relative to the 
magnitude cos a dz/dg and similarly along the characteristic 
tj = constant, the relation (2.15) may be considered as an ordi- 
nary differential equation relative to cos a ds/drj. By setting 
z = cos a dz/dg, the following 3-inear equation of the first order, 
based on equation (2.14), is obtained for determining z along 
the characteristic | = constant: 


dz + 1-Kx 1 ( a) z _ 1+a 1 ( g) dz _ 0 
dq sin 2a "2 sin a drj 


(2.16) 


The solution of this differential equation has the form 


z = exp 



1+q* (a) 

sin 2a 


L 





1 +a* (a) 

sin 2a 


*1 


1+a* ( g) ds 
2 sin a d£ 


(2.17) 


where tJq is the value of tj at a certain initial point on the 
characteristic considered, g = constant, and z Q = cos OQ(dz/dg) 0 

is the value of z at this point. The integral under the sign of 
the exponential function in equation (2.17) can be readily computed. 
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if the magnitudes entering the expression under the integral are 
expressed in terms of the Mach angle a, by making use of the 
relation 


X 



1 

J sin 2 a + (x-l)/2 


and the relations 

dt} = da = - d0 ^ = X tg a 

vhich hold along the characteristic of the first family £ = constant, 
the first of which follows from relations (2.7) and (2.11) and the 
second of which was already obtained in section I, part 4. The 
following expression is thus obtained: 


T] 



1 + a'(q) 
sin 2a 


dn = In 



sin 2 a [sin 2 clq + l/2(x-l)] 


sin 2 ctQ [sin 2 a + l/2(x-l)| 

J 

(2 


Hence, by replacing in equation (2.16) the magnitude z by its 
expression, and denoting by G(a,OQ) the expression under the 

logarithm sign in (2.18) 


cos a 


d l 

n 


1 

G(a,ao) 




1-hx 1 (a) ds 
2 sin a dt] 


drj 

(2.19) 


An entirely similar solution of equation (2.15) gives 


cos a 


ds 

dri 


1 



£ 

r .. X l+a*(a) dZ 
+ J G(a fao ) g -J&t df 

&0 


( 2 . 20 ) 


By making use of the property of mono tonicity of change of the 
magnitude of the velocity and the angle of inclination of the 
velocity vector along the characteristics in the local supersonic 
region (sec. I, pt. 3), it is noted that d2/d£ ^ 0 and ds/di] > 0. 
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Moreover, a 1 (a) > 0. 
the following inequalities 


Hence, from equations 
are obtained: 


(2.19) and (2. 20), 


cos a 



cos 



1 

G(a,a Q ) 


( 2 . 21 ) 


cos a 


ds 

di} 


cos 



1 

G(a,a 0 ) 


( 2 . 22 ) 


In order to obtain the accelerations at a certain point K of 
the local supersonic region, the characteristics of the first and 
second families are drawn from this point up to the intersection 
with the contour at points A and B, respectively (fig. 25). 

Let a = and a = a^g at points A and B, respectively. 

Also, let the radii of curvature at A and B be B^ and B g , 

respectively. 


The inequalities (2.21) and (2.22) at the point K are con- 
sidered where point B is taken for the initial point in inequal- 
ity (2.21) and in inequality (2.22) point A is used. By making 
use of relations (1.14), (2.1), and (2.7), approximations are 
easily obtained for the magnitudes d £/dZ at point B and drj/ds 
at point A 



B 



cos <Xq 2 


ds 



COS 


(2.23) 


By taking the direction of the x axis at point K along a 
streamline, the following equations are obtained: 


de 1 /d| + drA 

chE * 2 cos a \dZ ds/ 


^ _ .Xtg a /d T) _ d|\ 

2 cos a \ds d ij 


|£ i (si . ei\ 

oy 2 sin a \dZ d sj 


|X = _ S'- (SL i + * 1 ) 

oy 2 sin a \dZ ds/ 


(2.24) 

where the relation a* ( X) = - ctg a/ X is used. By making use of 
inequalities (2.21), (2.22), and (2.23) and equations (2.24) at 
point K, the following inequalities holding at this point are 
finally obtained: 
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H | S ^“oi'V + L ^ a 02 #E 2^] a 

If | - f(«01> E l) + »(«)-/ ctg a 

if I - + l Ccc 02 ,E 2 5] N(a) X tg a Vtg a 

Sxl r~ 

5y| - J^( a Ol> E i) + N(a)X*/tg a (2.25) 

where 


N(a) 


sin 2 a 


-|2(X-1) 


[sin <x + l/2(x-l)J 


L CctQ^E) = ~ /ctg cfo 


T £® inequalities (2.25) show that the absolute values of the 
derivatives cannot exceed certain finite values, depending on the 
Mach number at the given point of the supersonic region a nd the 
Blniaum filed radius of currature on the part of thfconSi^at 
is a boundary of the supersonic region. 

, a ,.^ he assumption that the cause of the formation of a shock wave 

if ° f ^ charact eristics within the local supersonic 
region is therefore invalid. Thus, it is possible to fo rmula te 

of 8 * Z ~: f0r c ertain conditions determining the mot ion 

of - the gas a poten tia l flow existed wit h loca l ?i£i55 

— fonaat ion of a shock wave within or on the boundary of the 

i ^personic re gion, arisi ^from any change in 

-i r! n i!i. | PreC9d ed by tne occu ^ rence of infinite accelerati55i~a t 
~ h interior point s of the supersonic region. Th at is . Tf~hh£ 

ln the gupersonic^region 


^When the flow about 
determining the motion of 
ing flow and the shape of 


the body is infinite, the conditions 
the gas are the Mach number of the oncom- 
contour of the body. 
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The theorem obtained refers only to the internal points of the 
supersonic region. From inequalities (2.25), no conclusion can be 
drawn as to the finiteness of the accelerations on the transition 
line because the right-hand sides of three of them approach infin- 
ity as the transition line is approached. Hence, the question 
whether the shock waves are preceded by the occurrence of infinite 
accelerations at the points of the transition line remains open. 


3. Criterion of Collapse of Potential Flow 

In section II, part 1, it was shown that a deformation of the 
contour can lead to the collapse of the potential flow. Before the 
deformation of the curvilinear segment of the contour at its points, 
the inequality &X/d(-e^.) > - Xtg a. holds. At the end of the 

deformation, however, at the points of the obtained straight seg- 
ment, dV^-^) = - oo, as follows from theorem 6. 

Hence, for ary intermediate state of the deformed contour 
there is first attained, at some point of the segment, the equation 

= - Xtg a (2.26) 


This state of the deformed contour is, in a certain sense, 
critical because for further deformation of the contour, inequal- 
ity (1.13) breaks down and thus the flow becomes impossible with- 
out a shock wave. 

In the given case condition, equation (2.26) is the criterion 
of breakdown of the potential flow. Equation (2.26) may also be a 
criterion for the breakdown of the potential flow in the case where 
the flow is about a fixed contour but where the Mach number of -the 
oncoming flow increases. 

If for a certain Mach number, a potential flow existed with 
local supersonic region and at the points of the contour, the 
inequality 


dX 


> 


X tg a 


were satisfied, then if with increase in the Mach number of the 
oncoming flow inequality (1.13) starting from a certain Mach number 
breaks down, this limiting Mach number would correspond to the 
attaining of equation (2.26) at some point of the contour. 
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The Mach number of the oncoming flow, for which at seme point 
of the contour condition (2.26) is satisfied, will be denoted as 
the breakdown Mach number M* . 

The number M* is thus the limiting Mach number beyond which 
the body is subject to wave resistance. 

The breakdown criterion (2.26) can be conveniently repre- 
sented in another form, by considering the magnitude 

6 1 * = 0 k “ Cp( \ ) (2 * 27) 


where cp(A) is a known function shown in figure 4. By differen- 
tiating this equation along the contour 


ds 


1 




fjfk 

ds 


But cp'(A k ) = l/(A k tg ct^) 


so that the following equation for 


the breakdown criterion is obtained equivalent to (2.26): 



(2.28) 


It should be noted that the breakdown of inequality (1.13) at 
any point of the contour does not lead to the breakdown of the 
potential character of the flow near this point but' makes impossible 
the ending of the characteristics of the first family from the 
neighborhood of this point on the transition line as the law of 
monotonicity would otherwise fail to hold at the transition line. 

The characteristics of the first family must, in this case, merge 
in the shock wave. 

In order to determine the magnitude M* of a given contour, 
it is necessary to know how for an increase in the Mach number of 
the oncoming flow, the value of the velocity at the points of the 
contour changes for a potential flow with local supersonic region. 

At the present time, certain methods are known for the approxi- 
mate solution of the problem of the flow about a fixed contour for 
Mach numbers exceeding M^. (See, for example, references 4 and 5.) 
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In all these methods, the convergence of the process of the 
successive approximations has not been proven j hence, the results 
of the computations made with the aid of these methods must be 
regarded with reserve. It seems probable, however, that up to 
the time when the potential flow is actually possible, these methods 
at least qualitatively represent the true character of the change 
in the distribution of the velocity on a fixed contour with increase 
In the Mach number of the oncoming flow. 

It is thus of interest to see whether from these solutions 
there is a tendency toward the attainment on the contour of the 
condition (2.28) with increase in Mach number. 

From the geometric data of a wall and distribution of the 
velocity along the wall, obtained in the work of Gortler (refer- 
ence 5) for a Mach number of the oncoming flow M = 0.9, the angles 
of inclination were computed using equation (2.27). The 

pattern of the lines of flow and the lines of equal velocities for 
this case is shown in figure 26. The dependence of 0-j_* on the 

coordinate along the wall % is represented in figure 27. 

As may be seen, the flow is near the breakdown. The monotonic 
decrease of the angle proceeds up to f » 5.5, after which the 
angle remains almost constant or slightly decreasing. 

In figure 28 Is shown the dependence of 0^ on x for the 

flow about the profile, which Is considered in an American Inves- 
tigation as M = 0.75 and M = 0.83. The velocity distribution 
over the contour and the boundaries of the local supersonic region 
are shewn In figure 29. In figure 28, it is very clearly seen 
that the tendency toward a breakdown of the monotonicity on the 
transition line increases with increase in the velocity. At 
M = 0.75, somewhat exceeding M cr , the angles decrease almost 
as fast as on the contour. For M = 0.83, the monotonicity of 
change of the angle 0-^ has already broken down. Hence, at 

M < 0.83 there is already a breakdown of the potential flow near 
this profile. Strictly speaking, the breakdown for M = 0.83 no 
longer has significance. Thus, in the given solutions, a tendency 
is revealed to attain on the contour condition (2.28) on increas- 
ing the Mach number of the oncoming flow. 
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4. Infinite Accelerations in Accurate Solutions 

In section II, part 2, it is shown that an infinite accelera- 
tion on a contour can occur only if the curvature of the contour 
at any point is infinite. 

An exact solution of a certain flow with local supersonic region 
with the occurrence of infinite acceleration at two symmetric points 
of the contour was given hy Ringleb (reference 6). According to the 
results of this section, this flow indicates an infinite curvature 
of the contour at these points. Figure 30 shows the flow pattern 
and the velocity distribution. At points P and Q of the contour, 
Ringleb obtained infinite values of the derivative of the velocity, 
that is, infinite accelerations. It is found that in the plane of 
the hodograph, at the points corresponding to P and Q in the 
flow plane, the transform of the streamline and the epicycloid have 
a common tangent. 

In the work of Karman during 1941 (reference 3), the Ringleb 
solution is analyzed in detail with the object of showing for this 
type of flow the possible reasons for the formation of shock waves 
in the local supersonic regions. Karmen showed that if at the con- 
tour in the supersonic region there is a point with infinite accel- 
eration and a value of f <» (where \Jr Is the stream func- 

tion and 0 is the angle of inclination of the velocity vector) 
then, as in the case of Ringleb, the transform of the streamline 
in the plane of the hodograph touches the epicycloid at a point 
corresponding to the point with infinite acceleration in the flow 
plane. K&rman identifies this contact in the hodograph plane with 
the presence of an infinite acceleration at the corresponding point 
of the flow plane. He does not, however, observe the fact previously 
pointed out that the presence of an infinite acceleration at any 
point of the supersonic flow indicates the presence In it of an 
infinite curvature of the streamline, and associates the infinite 
acceleration, In this particular case, with the impossibility of a 
continuous flow In the general case. This lack of observation is 
indicated by the fact that with the aid of the previously mentioned 
condition of the tangency in the plane of the hodograph, Kerman 
seeks an infinite acceleration on the contour of the NACA-4412 
profile, which has no infinite curvatures, making use of the exper- 
imental distribution of the velocity on the contour. In the plane 
of the hodograph, Karman obtains the point of tangency of the 
transform of the contour with the epicycloid but corresponding to 
this point of the contour of the profile, infinite acceleration 
does not occur as should be the case. No satisfactory explanation 
of this fact contradicting his results is given. 
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This noncorrespondence is explained by the fact that the tan- 
gency in the plane of the hodograph can occur simultaneously with 
the condition cty/dS = » at the point considered on the flow plane. 
An example of this is the Meyer flow, that is, the flow about a 
certain contour for which the entire contour coincides in the plane 
of the hodograph with the epicycloid but nevertheless, at the points 
of this contour, the acceleration is finite if the curvature is 
finite. 

Moreover, under the conditions of the problem of the flow with 
local supersonic region, the tangency in the plane of the hodograph 
of the transform of the contour with the epicycloid is equivalent 
to condition (2.26) at the point of the contour corresponding to 
the point of tangency, because at this point of the contour the 
acceleration is finite if the curvature of the contour is finite. 

In recent times, investigations have appeared on exact solu- 
tions of the problem of flow with local supersonic regions where 
the equations of S. A. Chaplygin, in the plane of the hodograph, 
are used. Notwithstanding the fact that in these papers examples 
are given of the existence of potential flows with local supersonic 
region, the value of these investigations is limited by the fact 
that with change in the Mach number of the oncoming flow there 
occurs simultaneously a deformation of the contour of the body. 

Due to this fact, points with infinite acceleration appear within 
the supersonic region, which on the basis of the theorem proved in 
section II, part 2 indicates an infinite curvature of the contour 
at these points. 

From the foregoing considerations, it Is evident that it is a 
mistake to associate the occurrence of such a limiting line of the 
contour with the impossibility of a potential flow with local 
supersonic region about a fixed profile. 

The Investigation of the principal properties of a flow with 
local supersonic region must, in the opinion of the author, be 
carried out directly in the flow plane. 

In the work described herein, all results were obtained 
directly from investigations in the flow plane. 


Translated by S. Eeiss 
National Advisory Committee 
for Aeronautics 
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Figure 19. 
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Figure 29. 
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Figure 30. 



